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Abstract 

In this paper, the effects of a newly presented structure of 
unified power flow controller (UPFC) which is based on two 
shunt voltage-source converters and a series capacitor, are 
investigated on the dynamic stability of power system. A 
power system including the newly presented structure of 
UPFC is introduced to examine this. The non-linear dynamic 
model of the system under investigation and its non-linear 
equations are presented as well. The linear equations are 
achieved by linearizing the non-linear equations around the 
operation point. A modified Heffron-Philips model is 
proposed through the linearized equations of the 
investigated system. The state equations of the system are 
then extracted using the modified Heffron-Philips model. 
The damping controller of UPFC is then presented and its 
control parameters are calculated by phase compensation 
method. Finally, the most effective input signal of UPFC 
control signals is selected using some techniques. The 
validity of the theoretical investigations is confirmed by 
simulation results obtained in MATLAB/Simulink. 

Keywords 

UPFC Modeling; FACTS Device; Electromechanical Oscillations 
Damping; Low Frequency Oscillations; Modified Heffron-Philips 
Model 

Nomenclature 



M 


Inertia constant (second) 


X Eg for g = 1,2 


Transformer g reactance (pu) 


X -el for k = a,b,c 


Equivalent reactance of capacitor in k phase 

(PU) 


X t for i = d,q 


Steady state synchronous reactance of the 
generator on i axis (pu) 


X 4 


Direct axis transient synchronous reactance of 
the generator (pu) 




Reactance of transmission line (pu) 


X le 


Reactance of transformer (pu) 


c, 


Dc link capacitor (pu) 


c. 


Series capacitor (pu) 


R E t for g = 1,2 


Equivalent resistance of shunt converter g (pu) 


I b for i = d,q 


Current in infinite bus on i axis (pu) 



I tl for i~d,q 


f^nrrpnt nf ntiaQP k in fr^TTiTniQQinn linp nn i 

oAlo I Lf LI 1 


I ut for i = d,q 


V_QILcLll KJl UlldbC /C 111 SCLICS CdUdCLlUI OIL ( dAlO 
\P U ) 


1 for o-l 2 

i = d ,q 


Current in converter g on axis (pu) 




Voltage at dc link (pu) 


V lel fork = a,b,c 


Series capacitor voltage of k phase (pu) 


v b 


Voltage at infinite bus (pu) 


V, 


Voltage at generator terminal (pu) 


V tt fori~d,q 


Terminal voltage on axis (pu) 




Reference voltage (pu) 


AW 


Provided electrical torque by UPFC damping 
controller (pu) 


P 


Input mechanical power to the generator (pu) 


p 
/ 


Electrical power of the generator (pu) 


A 


State matrix 


B 


Input matrix 


C 


Output matrix 


D 


Direct transfer matrix 


X 


State vector 


Y 


Output vector 


(J 


Trrnni - vprfnr 

J_l LLJ U. L V CvAVJl 


F(s) 


Forward path of the UPFC damping controller 


K c (s),K (s), 
K L (s),KJs) 


UPFC damping controller transformer controls 


H(s) 


Transfer function of UPFC damping controller 


/, 


Switching frequency (Hz) 


S 


Phase angle of System load (degree) 




Phase angle deviation of System load (degree) 


S E, g for g =1,2 


Phase angle of shunt converter g (degree) 


g for g = 1,2 


Phase angle deviation of shunt voltage source 
converter g (degree) 


m E, t f or g =1.2 


Pulse width modulation index of shunt 
converter g 


Am E .« M g=l,2 


Pulse width modulation index deviation of 
shunt voltage source converter g 


« 


Fundamental angle frequency ( rad s~' ) 




Oscillation angle frequency ( rad s~ l ) 


Aco 


Error signal ( rad s~' ) 


D m 


Friction factor of generator 
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^UPFC 


Friction factor of the UPFC damping controller 




Input signal of the UPFC 




System operation domain 




UPFC controller gain 


K a 


Gain of AVR in generator 


K t - K a 


Constants of Hef fron-Phillips model 


K cui ,K pui ,K iut , 


Constants control of input signals 


K , , K . , K , 

qd ' pi! ' vd 


Constant control of input signal 


r 

Jo 


Direct axis open circuit time constant of the 
generator (second) 


T 


Time constant of the generator of AVR 




Time constants of phase compensator 



Introduction 

Nowadays, considering the increasing growth of 
power systems in different countries and their mutual 
connections and also the great interest in transmitting 
the possible maximum power through the transmission 
lines, the maintenance and enhancement of power 
systems stability becomes more important [Shahir 
(2013), Kundor (1994)]. The modern power systems are 
designed to achieve high generation capacities, 
bringing more risk into power systems stability. One 
of the most important factors threatening power 
systems stability is the electromechanical oscillation of 
the generator. Dynamic instability is raised due to 
balance lose between the output electric power and 
input mechanical power, lack of damping torque, and 
sudden load distortion outbreak. The frequency of 
such oscillations falls in the range of one to few Hertz 
[Shahir (2013)]. Several factors such as operation 
condition, load characteristics, lines impedance and 
power, series capacitors, and voltage regulators 
influence this type of oscillations. 

One of the most effective solutions is to install power 
system stabilizer (PSS) in generator location to damp 
the electromechanical oscillations [Yassami (2010)]. 
The PSS generates a complementary signal to rapidly 
damp the distortions caused by oscillations of power 
system. The PSS usually considers frequency error 
signal as its input and generates a signal added to the 
voltage. The practical and simulation results show that 
the PSS will be highly effective in system oscillations 
damping if it is properly regulated. 

Recent technology development in semiconductor 
based devices and their control techniques field and 
high interest of power systems designers in utilizing 
flexible ac transmission systems (FACTS) for maximum 
power transmission has provided a new stability 
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research field. The application of Applying FACTS 
devices brings about more reliability and less costs in 
different network arrangements under distortions 
occurrence conditions [Song (2004), Xia (2010)]. Among 
the FACTS devices, the voltage-source based converters 
considerably affect the transient state and improve 
dynamic stability of power systems [Hingorani (1999)]. 
The UPFC, one of the most important and effective 
FACTS device applied in series-shunt mode consists of 
a series and a shunt three-phase two-level converter 
connected to each other through a dc capacitor. The 
series converter controls the active and reactive powers 
of transmission line by controlling the amplitude and 
phase angle of series voltage and the shunt converter 
compensates the UPFC losses and stabilizes dc link 
voltage in addition to transmission line reactive power 
control [Sarvana (2009)]. Many different structure have 
been proposed for UPFC [Shahir (2012), Sadigh (2010)]. 
In [Sadigh (2010)], a new UPFC structure consisting of 
two shunt voltage-source converters and a dc capacitor 
installed between the mentioned converters is 
presented in which the UPFC parameters are controlled 
by the series capacitor. The voltage and the current of 
the series capacitor are controlled by the operation of 
the two shunt converters [Sadigh (2010)]. The control 
and the measuring systems of conventional UPFCs are 
separately designed due to the different insulation and 
protection processes of series and shunt converters. This 
leads to high time and cost expenditure of such devices 
design and manufacturing. In the structure of UPFC 
presented in [Sadigh (2010)], the evaluation bulk is 
decreased due to similar shunt converters application 
and as a result; the control and measuring systems 
design and manufacturing costs are significantly 
decreased. Also, in comparison with congenital UPFC, 
total VA rating for high active and reactive power flow 
control is decreased. In the structure of [Sadigh (2010)], 
the signals could be applied as UPFC input signals for 
stability and electromechanical oscillations control that 
is an advantage of that structure. Three signals can be 
applied for stability control in conventional UPFCs 
[Shahir (2012)], while four signals are used to control 
the operation of the newly presented UPFC. Reduction 
of the THD value of the injected series voltage and 
applying series capacitor as constant series capacitive 
compensator during converter 2 standby interval are 
other advantages of the structure presented in [Sadigh 
(2010)]. 

The modeling process of conventional UPFCs and its 
effect on dynamic stability is investigated in [Fang 
(2003)] and transient stability rate improvement is 
evaluated in [Sreenivasachar (2000)]. In this paper, the 
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effects of the UPFC structure presented in [Sadigh 
(2010)] on the dynamic stability and electromechanical 
oscillations damping of power system are investigated. 
The analysis starts by nonlinear equations extraction. 
The linear equations and the modified Heffron-Philips 
model are obtained by linearizing the nonlinear equations 
around the system operation point. The phase 
compensator technique [Larsen (1995), Chandrakar 
(2010)] is then applied and the proper controller is 
designed and regulated. Finally, the most effective 
input signal for UPFC control is obtained by simulating 
the system in MATLAB/Simulink to damp the 
oscillations of synchronous generator's rotor. 

Investigated System 

A sample power system including a single machine 
with infinite bus along with the UPFC presented in 
[Sadigh (2010)] is used to investigate the dynamic 
stability. The single line structure of the investigated 
power system is illustrated in Fig. 1. The utilized 
UPFC consists of two parallel transformers, a series 
capacitor and two shunt voltage-source converters 
controlled by pulse width modulation (PWM) technique 
[Shahir (2012)]. Here, R El and R E2 are neglected in 

comparison with X E l and X E2 , respectively. 
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FIG. 1 THE SINGLE LINE DIAGRAM OF INVESTIGATED SYSTEM 

Dynamic Modeling of the Investigated Power 
System 

Nonlinear Dynamic Model 

All reactance and resistances of the generator, 
transformers, and transmission line should be 
considered for nonlinear dynamic modeling of power 
system. Therefore, the following nonlinear equations 
are obtained assessing the system shown in Fig. 1: 

8 = co Q (co-l) (1) 



co ■ 



P m -P e -D„Aco 



F' =- 



M 

-E„+E 



fd 



(2) 
(3) 



-fd 



-E fd +K a (V ref -V t ) 



3m 



E.\ 



* dc ■ 



AC, 



dc 



[cos § E \ sin 8 E1 J 



3m E 2 r 

H ~\ cos ^e,2 s i R ^E,7 

4C, 



dc 



l E\,d 
^E,l,q 
T E,2J 
l EXq 



V se =—I seJ M i = d,q 

The folio wings are valid for (l)-(6): 

P e =V t4 I t4+ V t , q I t , q ,E q =E' q +{X d -X' d )I t4 ; 

v, =y t ,d + jv t , q ;v t , d = x q h, q y,, q = E' q -X' d I,j ; 

h,d = he,d + l E\,d ! h,q = he,q + 1 E,\,q > 



(4) 



(5) 



(6) 



he,d =a x V b cosS- 



-cos S E 2 + c x E' q - 



cos 8 El 



W ■ X , b l m E,2 V dc . r , d \ m E,^dc . r 

he, q = a Vb sm 8 + sin 8 E _ 2 + sin 8 EA 



hxd =a 2 V b cosd + - 



b 2 m E,2 V dc 



cos<5 £2 +C 2 E' q + 



-cosS El 



hxq = a 2 V b sm S + - sin S E2 + sin 8 E1 

,/ • X b l m E,2 V dc s d l m E,\Vdc x 

h,2,d = a 3 V b sm 8 + cos<y £ 2 + c 3 £ + cos 8 EA 



'E,2,q ~ u 3 y b 



a-,V h sin 8 + - 



h m E,2 V dc 



sin 8 E2 + 



d l m E,\Vc 



dc 



sin 8, 



EA 



m E 2 V dc 

I t d =(a l + a 2 Wb cos <5 + (b\ +b 2 ) — '- cosS E 2 

l \r> /J.J , m E,l^dc o 

+ (q +c 2 )E q + (d t + d 2 ) — - — cosS El 



m E,2 V dc 



sin 8 



E,2 



l t, q =(a l +a 2 )V h s\x\S + (b x +b 2 )- 

m F \^dr 

+ (d x +d 2 ) E ^ dc sin<? £t i 



XbB = %E,2 + Xbv'Xal = X.w + ^E,l^te + ^te^se' 

X a2 = X BB X E y,X a3 = X E2 X El ;X hl = X E1 +X te ; 

= X E,1 +X te> X b2 = X Et 2X bv ;X b3 = X BB X se ; 



^i - X al X a2 + X bl X b2 ;X 2 - X bv + 



X p i 



The other constants are presented in Appendix A. 
Linear Dynamic Model 

The linear dynamic model of the sample system is 
obtained linearizing (l)-(6) around the system operation 
point. The following equations are obtained after 
linearization: 

A8 = co Aco (7) 
AP m -AP e -DAco 
M 



Acb ■ 



106 



Electrical Engineering Research Vol. 1 Iss. 4, October 2013 



www.seipub.org/eer 



-AE a +AE fd 
AE' = q 1 



T' 

l do 



AE fd =- 



-AE fd+ K a (AV ref -AV t ) 



AV dc = K,A8 + K,AE' q - K 9 AV dc + K cEX Am EX 

+ K c,5E^ 5 E,\ + K c,E,2 Am E,2 + K c,5E,2 Al 
A ^ = K m A $ + K ll E 'g + K sl Am E,2 + K s2 Am E,l 



(9) 

(10) 
(11) 
(12) 

'(13) 
(14) 
(15) 



+ K n AV dc + K s3 AS E2 + K s4 AS EX 
The followings are valid in (7)-(12): 

AP e = K,AS + K 2 AE' q + K pEl Am EX + K pSEl AS EA 

+ K p,E,2 Am E,2 + K p,SE,2 AS E,2 + K pd A ^dc 

AE q = K A A8 + K 3 AE' q + K qEl Am EX + K qSEl AS E1 
+ K qtEa Am E2 + K qSE2 A8 E2 + K qd AV dc 

AV t = K 5 AS + K 6 AE' q + K vEX Am EX + K vSEl AS EA 

+ K V ,E,2 Am E,2 + K v,SE,2 A5 E,2 + K vd A ^dc 

The constants used in (11)-(15) are presented in details 
in Appendix B. 

The Proposed Modified Heffron-Philips Model 

Here, it is possible to propose the modified Heffron- 
Philips model of the system shown in Fig. 1 based on 
(7)-(12) to analyse and investigate the system stability 
against oscillations and to evaluate its small-signal 
stability [Shahir (2012), Kundor (1994)]. The proposed 
modified Heffron-Philips model is a linear model 
applied in low frequency oscillations and small signal 
stability investigations. Fig. 2 shows the block diagram 
of the proposed modified Heffron-Philips model 
resulted from (7)-(12). In this model, one of the 
linearized control signals of the sample power system 
including UPFC, Am El , Am E2 , A8 E , , and A5 E2 is 

applied to the UPFC and controls the investigated 
UPFC's operation. 

The block diagram of proposed modified Heffron- 
Philips model shown in Fig. 2 contains 35 constants 
that can be defined as row vectors of K uk , K q uk , 

K v,uk , K c,uk > and K se ,uk as follows to simplify the 
calculations: 





= ^ K p,E,\ 




Kp,E,2 


K P,SE,2\ 


(16) 


K q ,uk 


= ^ K q,E,\ 


K q ,SE,l 


K q , E ,2 


Kq,8E,21 


(17) 


K v ,uk 


= \- K v,E,\ 


K V ,SE,\ 


K v,E,2 


K v,SE,2\ 


(18) 


K c , u k 


= l K c,E,l 


K C ,8E,\ 


K c,E,2 


Kc,SE,2i 


(19) 


K se ,uk 


= \- K s,E,\ 


K s,SE,\ 


K s,E,2 


K s,SE,2^ 


(20) 



The row vectors of (16)-(20), are applied on proposed 
modified Heffron-Philips system as the following 



input vector: 

U =[Am El A8 El Am E2 A8 E2 ] T 



(21) 




FIG. 2 THE PROPOSED MODIFIED HEFFRON-PHILIPS MODEL 
OF INVESTIGATED SYSTEM 

The State Equations of the Investigated System 

Based on (7)-(15), the dynamic model of the investigated 
system in the form of the state space equation can be 
expressed as follows: 

X=AX+BU (22) 
where X and U are defined as follows: 



X = [AS Aco AE' q AE fd 



AV dc AV S J 



U = [Am E ! A8 E 



Am 



E,2 



AS, 



E,2l 



(23) 
(24) 



A and B matrices of (22) are considered as follows: 
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K q,SE,2 






T d0 
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K a K v ,E,i 


K a K v,SE,\ 


K a K v,E,2 


K„K Vi se,2 
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T 
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T 
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K c,E,\ 




K c,SEA 


K c,E,2 


K Ct SE,2 


K s,E,l 




K s,SE,l 


K s,E,2 


K Sy SE,2 



(25) 



(26) 



It is possible to design an appropriate controller using 
(22) and select the most proper input signal of UPFC 
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from four input signals according to criteria detailed 
in the next section in order to maintain the system 
stability and to control the damping of power system 
against the rotor oscillations of investigated system's 
synchronous generator. Therefore, it is important to 
select and design appropriate input signal and 
desirable damping controller. 

UPFC Damping Controller 

The linearized model of the UPFC installed in power 
system explained in (22) can be shown in Fig. 3 [Wang 
(1997), Shahir (2012)]. The UPFC damping controller 
consists of two sections of input signal and UPFC 
controller. The input signal is applied on power 
system through the input transfer functions of K (s) , 
K L (s), and K c (s) . The state space equation obtained 
in (22) can be considered as follows [Wang (1997), 
Shahir (2012)]: 



~ AS 







C0 





~AS~ 









Acb 








~ a d23 


Aco 


+ 


~ b d2 


U (27) 






_ a d31 


a d32 


a d33 _ 


_ z m 




~ b d3. 





d\ 



~d2 



~d3 



AS 
Aco 

7 



+ DU 



(28) 



A3 



Electromechanical oscillation loop 
_ 1 Am 



K c {s) 



Au k 



FIG.3 THE CLOSED LOOP BLOCK DIAGRAM MODEL FOR UPFC 
DAMPING CONTROLLER 

In (27), k m and d m are the synchronizer and damper 
components, respectively, and z m consists of other 
state variables of (22). The oscillation angular 
frequency is obtained by co i = ^Jco Q k m . Here, the 

transfer functions presented in [Wang (1997), Shahir 
(2012)] are explained using (27) and (28). 

K c (s) is the transfer function from U input to Acb 
and obtained as follows: 

K c (s) = a d ^(sl-a d ^y x b^+b d2 (29) 

K (s) is the function of transferring rotor speed to the 
y measuring point shown in block diagram of Fig. 3 



and denoted as follows: 

Ko (*) = c*3 (si ~ a d33 )" 1 ^a d3l ^ + a d32 j + ^yC dl + c d2 j (30) 

K L (s) is the function of transferring U input to the y 
measuring point shown in Fig. 3 and calculated as 
follows: 



K L( s )= c d3( sI - a d33) b d3+ D 



(31) 



K m (s) is the transfer function of AS to Acb and 
calculated as follows: 



K m(s) = K + — d m +a d23 (sl-a di3 ) 



-V31 



S 

+ — a 

CO, 



d32 
'0 J 



(32) 



H(s), as a UPFC controller transformer, controls the 
operation of the system shown in Fig. 1 for different 
input signal values. In this investigation, the UPFC 
damping controller is considered ideal. An ideal 
UPFC damping controller should provide a pure 
positive electric torque in oscillation loop, damp 
electromechanical oscillations, and stabilize the 
generator rotor oscillations. Considering Fig. 3, the 
provided electric torque by UPFC damping controller 
could be calculated as follows [Wang (1997), Shahir 
(2012)]: 

%c= M)W) Affl (33 ) 
l-K L (s)H(s) 

The torque magnitude provided by UPFC damping 
controller is expressed as follows: 

(34) 

oscillations. 



at; 



UPFC 



D UPFC Aco 



D, 



UPFC 



reduces the generator rotor 



Considering (33) and (34), D UPFC can be regarded as 
follows: 

K c (s)K (s)H(s) 



D, 



UPFC 



(35) 



(36) 



\-K L {s)H(s) 
It can be simplified as follows: 

Dupfc = (K c (s)K (s) + D UPFC K L (s))H(s) 
= F(s)H(s) 

F(s) in (36) varies with power system operation 
conditions and is defined as follows: 

F(s) = K c (s)K (s) + D UPFC K L (s) (37) 
The operation conditions of the power system, 
applying UPFC structure presented in [Sadigh (2010)] 
as a damping controller, depend on the generator 
operation point, power system condition, and the 
UPFC input signal. Considering this fact, the selection 
of UPFC input signals and the most proper of them 
should be in away that the system stability against 
electromechanical oscillations is bonded and it achieves 
its stability as soon as possible. Also, the system 
should maintain its stability as electromechanical 
oscillations occurwhile the damping speed in a power 
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system is highly important. As the damping time is 
long, the power system would get damaged. 

UPFC Damping Controller Design 

If it is assumed that a power system has the operation 
domain of , F(s) would be a function of 

operation condition of ju and UPFC input signal and 
the controller is designed and installed on power 
system on this base. The desired controller for 
oscillation damping is the one that can guarantee the 
system stability, can acceptably react against the 
generator rotor oscillations, and can damp such 
oscillations in the shortest time. Proper input signal 
selection is one of the most important steps in the 
design of an appropriate controller. For a good 
controller designing, several approaches have been 
evaluated. In this paper, the phase compensation 
method is applied [Chandrakar (2010)]. For proper 
control signal selection of UPFC damping controller, 
three techniques presented in [Wang (1997), Shahir 
(2012)] are applied in the following subsections. 

1) First Technique 

The parameter based proper control signal 

selection for UPFC damping controller design is as 
follows: 

^selected = min ^ F ( s > M e ^(/O ( 38 ) 

Here, the selection aims to making it possible to 
design UPFC damping controller in a way that the 
system would be able to maintain its stability 
against rotor oscillations under any circumstance. 
This is accomplished targeting proper control 
signal selection in order to obtain the power system 
stability in minimum value of . 

2) Second Technique 

One of the most important issues considered in the 
designing process of damping controllers is to 
achieve the maximum control domain with 
minimum expense. Therefore, the proper input 
signal is the one that has the maximum effect on 
UPFC damping controller. Proper input signal for 
UPFC damping controller design on the basis of 
achieving minimum control system costs and 
maximum control signal effect is shown in (39): 

U selected = maX u t . ^ ? ( S ' M 'selected i u k) 



11 K £ m E,l' m E,2'^E,l'^E,2 



(39) 



3) Third Technique 



operation domain, ¥(//) . This is done under 
several conditions of to guarantee the power 

system condition and tries to select an input signal 
that would not considerably vary as varies 
and has a constant and invariable reaction against 
generator rotor oscillations in ¥(//) domain. If the 
influence of damping is increased, other system 
parameters and its operation condition under 
different values might be undesirably 

affected. On the other hand, significant decrease in 
the influence of damping controller can decrease 
system stability level and make the power system 
unreliable. In (40), a technique is introduced to 
achieve power system stability against rotor 
oscillations with the least impact on other power 
system parameters and components. However, in 
the cases, the results are not meaningful or not 
expectable; it is possible to apply (39) to guarantee 
the results obtained from (40) for input signal 
selection of damping controller designing. 

max^ F(s, /u,u k )- min /( F(s, /u, u k ) 



selected 



min F(s,M,u k ) 



(40) 



UPFC Controller Design 

As shown in Fig. 3, the UPFC damping controller 
consists of H(s) feedback transfer function. In [Shahir 
(2012)], the H(s) function is investigated in conventional 
UPFCs damping controllers applied in single-machine 
power systems. In this paper, by developing the 
method in [Shahir (2012)], H(s) is considered as 
follows: 



H(S)- 



sT u , 



K„ 



l+sT 2 l + sT 4 



(41) 



In third technique, (40) is applied to select the most 
appropriate input signal considering system 



\ + sT w l + sT u 1 + sTj l + sT 3 
Relation (41) can be expressed as follows: 

H(s) = \H\Z0 (42) 
where \H\ is the amplitude and 9 is the phase angle 
of H(s) . Considering (37), F(s) is rewritten as follows: 
F(s) = \f\Z</> (43) 
In (43), \F\ is the amplitude and <p is the phase angle 
of F(s) . Both H(s) and F(s) transfer functions 
considerably affect D UPFC and oscillation damping 
capability of power system (see (36)). 

In UPFC damping controller design, the controller 
tries to apply a pure positive torque on investigated 
system using the phase compensation technique. 
Based on (36), F(s) is a function of and it is 

possible to achieve an ideal controller that provides 
positive torque for generator oscillations damping, just 
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by controlling H(s) parameters. Considering (36) and 
(42), in order to provide a pure positive torque and to 
design on phase compensation technique, the 
following relation is valid: 



H 



D, 



UPFC 



Z0 = Z- 



(44) 



where the amplitude and the phase angle of H(s) is 
obtained and regulated by (44). 

Simulation Results 

One of the most important aims of UPFC installation is 
to control power system against low frequency 
oscillations happening as a distortion occurs in power 
system and the generator rotor is driven out of its 
synchronous mode. UPFC is in charge of proper 
responding as such oscillation occurs. Fig. 4 shows the 
power system illustrated in Fig. 1, where the Aa> error 
signal is sampled and applied to UPFC through U 
input signals. In the next parts, the proper signal is 
selected through analytic methods and UPFC 
controller regulation approaches are given. Nonlinear 
simulation is then accomplished on the investigated 
system. 



Aco 



V u + 
"X™ 



^aIw Excitation Excitation 
j j| Transformer 1 Transformer 2 \^ 




jYm_| 



UPFC 



Converterl Converter! 

An, 



damping 
controller 



FIG. 4 UPFC TO SIMB CONNECTION 

Proper Input Signal Selection 

F(s) is a function of UPFC damping controller input 
signals (37) and directly relates to and varies 

with and changes AT UPFC magnitude. Therefore, 

F(s) can be applied as a parameter of UPFC controller 
signals evaluation and the most proper control signal 
recognition. The most proper UPFC controller signal is 
selected utilizing techniques introduced in previous 
section and investigating F(s) under different 

conditions considered as Table 1 in this paper. The 
magnitudes and the angles shown in Table 2 are 
obtained evaluating F(s) under jj. x and ju 2 operation 
condition applying Nichols [Wang (1997), Shahir 
(2012)] approach. Fig. 5 shows the oscillation 



frequency in ^(ju) . It is initially observed that the 
control system is stable for all *¥(ju) values if the control 
system is designed under /j 2 operation condition. 
However, this requires further attention especially 
about the UPFC device applied in this paper whose 
operation is controlled by four input signals of Am E , , 

Am E 2 , AS E j , and A5 E2 . In conventional UPFCs applied 

in power systems, the number of input signals is 
decreased to three signals and the investigations are 
carried out on these three inputs due to the incapability 
of using A8 b [Wang (1997)], the phase angle deviation 

rate of the series voltage-source voltage. 

TABLE 1 SYSTEM OPERATION CONDITION IN T(//) DOMAIN 



operation condition 


V : 


v b 


P e 


M, 


l.Opu 


l.Opu 


OXpu 


M 2 


1.0 pu 


l.Opu 


l.Opu 



TABLE 2 THE AMPLITUDE AND PHASE ANGLE OF F(s) UNDER Mi AND 
jU 2 OPERATIN CONDITION 



«J 


\F(s)\ toclB 


phase angle to degree 


Mi 


Mi 


M 


Mi 




32.3789 


419.6940° 


294.8861° 


0.2160 




134.9189 


430.3020' 


407.4034° 


12.1071 




14.8118 


124.6427° 


313.9723° 


1.4827 


S E,2 


41.7253 


507.4505° 


300.6187° 


0.3162 



FIG. 5 OSCILLATION FREQUENCY CREATED IN 
INVESTIGATED SYSTEM 

1) Case study 1 

In this section, the approach mentioned in (38) is 
applied to detect the controller signal. Based on 
(38), the signal that has the minimum value under 
/u x operation condition can be selected. Therefore, 
applying (38) and based on the results shown in 
Table 2 results for the most proper signals of 
controller design are m E l and S E 2 , respectively. 

2) Case study! 

Relation (39) selects the controller signal under n 2 
operation condition in a way that the input signal 
has the maximum amplitude. Therefore, based on 
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the results shown in Table 2, the proper signal for 
controller design is S E1 . 

3) Case study 3 

Considering (40) and Table 2, and selecting 

U k = m E,\ ' U k = $E,\ ' U k = m E,2 ' U k = 3e,2 ' aS 

different input signals, the optimal input signal 
("selected ) is obtained using u k =m El . 

UPFC Controller Parameters Regulation 

As mentioned before, the UPFC controller parameters 
are regulated applying phase compensation technique 
(relation (44)). The UPFC controller parameters 
displace the Eigen values and are regulated in terms of 
UPFC damping controller input signals in order to 
create desired D UPFC in for generator rotor 

oscillations damping and maintain system stability 
[Shahir (2012)]. As the system Eigen values are less in 
terms of a specific input signal in damping controller, 
the input signal is more effective in damping 
controller and as the UPFC controller gain is less, 
lower K v , the UPFC damping controller can be 
designed and utilized with lower cost. The UPFC 
controller parameters applied in this paper are as 
7] =0.9* , T } =0.9* , T u = 0.01s , and T w =Ws . Other 
UPFC controller parameters are obtained from (44) by 
input signal change. The system Eigen values and 
other UPFC damping controller parameters are 
respectively shown in Table 3 and Table 4 under n x 
and n 2 in terms of different input signals. Based on 
the results shown in Table 3, the most appropriate 
signals for UPFC controller input under jj x operation 

condition is 8 E x from effectiveness point of view and 

is m E 2 from the view point of cost. Based on the 

results obtained from system analysis in terms of each 
UPFC damping controller signals, the maximum 
damping is achieved selecting S El as the UPFC 

damping controller input signal while the minimum 
cost is achieved under /u 2 operation condition in 
which m E 2 is selecting. 



TABLE 3 UPFC CONTROLLER PARAMETER UNDER //, 



Input signal 


UPFC controller parameters 


Eigen values 


">Ei 


K u =0.9, T 2 =0.2, 7/ 4 =0.15 


s = -0.1302+ 70.0460 




K„ =2.3, T 2 =0.12, T 4 =0.10 


s = -0.4718 + ;0.3120 


m E2 


K u = 0.6, T 2 = 0.65, T A = 0.09 


s = -0.1368± yO.0101 




K u =0.9, T 2 =0.55, T 4 =0.20 


s = -0.2381± j'0.1167 
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TABLE 4 UPFC CONTROLLER PARAMETER UNDER fJ 2 



Input signal 


UPFC controller parameters 


Eigen values 


m EX 


K u =\.\, T 2 =0.95, r 4 =0.95 


s = -0.0470 + y'0.0747 


s E , 


Jsr u =o.65, r 2 =o.25, r 4 = o.i5 


s = -0.2074 + 7'0.3264 


m E.l 


£,,=0.3, r 2 =0.72, J 4 =0.65 


s = -0.0727 + ;0.0241 


4.2 


K u =3.3, T 2 =0.65, T 4 =0.45 


s = -0.0903+ y'0.0003 



Nonlinear Analysis 

In this section, the system shown in Fig. 1 is simulated 
under relations mentioned in (l)-(6). The aim is to 
investigate the way system responds to sudden load 
variation caused by distortions. The simulation is 
carried out under the circumstance that a three-phase 
to ground fault occurs in the system shown in Fig. 1. 
The fault occurs in t = 1 sec . In this simulation, the 
operation condition of the system is as Table 1. Other 
simulation parameters are detailed in Appendix C. Fig. 
6(a) shows the simulation results under operation 
condition. The results are evaluated for all four input 
signals. The system condition during distortion occurrence 
without applying UPFC damping controller is shown 
in Fig. 6(a). 



14 




45 1 1 1 1 

2 4 6 8 10 

(b) 

FIG. 6 SIMULATION RESULTS (a) UNDER //, (b) UNDER fj 2 

Under such condition, the system starts to oscillate 
and the stability is weakened. Also, Fig. 6(a) shows the 
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system condition under distortion occurrence and at 
the presence of UPFC damping controller designed for 
m E j , 8 E j , m E 2 and <J £ 2 input signals, respectively. 

According to the results shown in the figure, the 
system stability is reinforced and the oscillations are 
rapidly damped under /u x and applying u k = 8 El , 

while the system damping is slow under u k = m E 2 

condition, which leads to the defection of the power 
system equipments. Fig. 6(b) shows the system 
response to the created distortion when it operates 
under ju 2 operation condition. In this figure, it is 
illustrated that the investigated system oscillates if no 
UPFC damping controller is applied. Also, Fig. 6(b), 
shows the assessment results of system damping rate 
improve the application of UPFC input signals and the 
damping controller for m El , 8 El , m E2 and 8 E 2 , 

respectively. According to the results shown in Fig. 
6(b), the system is more stable against the oscillations 
with u k = 8 E1 . Also, it is not reliable with u k = S E 2 

and the oscillations are damped very slowly. 
Conclusions 

Although several approaches are used to select proper 
input signals, different standards can be considered to 
select such signals for UPFC damping controller 
installed in power systems (such as the system 
investigated in this paper). If the purpose of damping 
controller design is just to control generator rotor 
oscillations and the economical concerns are not under 
consideration, the most appropriate UPFC damping 
controller signal is the one applied to system desirably 
in all and damps oscillations rapidly and 

controls the system. Therefore, it is possible to utilize 
m E 2 as input signal. If the economical concerns and 

costs reduction interests are considered in power 
systems designing process, the proper controller is the 
one that possesses lower UPFC controller gain for 
specific input signal in all ^(ju) domains. Therefore, 
S E j can be the most appropriate option. The stability 

of the controllers designed in power system is 
significantly important. The stability of the designed 
controller can be explained as such statement that the 
damping controller input signal is effective in all 
domains and operates invariably. Considering 
such a standard, selecting m El can guarantee the 
designed controller's stability in all domains. 

APPENDIX A 

The constants used in (l)-(6) is as follows: 



X a3^- X te ~ X E,0 



x \ x e,\ 



1 



(l-jX 2 X a3 ); b x = 



r -x x } 

X l X E,l , 
( X E,1 ~ X te)( X a3 + X al) 



J X BB 

~X te (X a3 -X a2 ) 



x \ x e,\ 



X,. 



X,X, 



f X ^ 
V X E,i J 



(X a2 + X b3 + X b2 )(X Erl -X, e )+ jX, 



c, = 



X 1 X E,1 



x 



bv 



X, 



jX 2 (X a3 +X b3 +X a2 ) 



(X te -X E , l )(X h3 + X b2 )-jX l 

X l X E,l 

~-X te (X b3 + X b2 ) + jX 1 



d,=- 



X l X E,l 



-jx 2 (x b3 + X b2 )- 



f x h ^ 



X, 



APPENDIX B 

The constants used in (11)-(15) are as follows: 

K l =[(a l +a 2 )(V t ^+X q I td W b co & S 



-[(a l +a 2 )(V tJ -X' d I tq )]V bS mS 
K 2 = (V t , d -X' d I t , q )( Cl +c 2 ), K 6 =(c l+ c 2 ) + 



K 3 =l + (X d -X' d )(c l +c 2 ), K 4 =-(a 1 +fl 2 )y z ,sinJ 



t.d 



K 5 = (aj +a 2 )X' d V b sin J + (aj +a 2 ) 
3a 3 m E 2 V b 

K 7 = ' (cos 8 sin 8 E 2 - sin 8 cos 8 E 2 ) 

3a 2 m E x V h 

+ ■ — (cos 8 sin 8 E x - sin 8 cos 8 E x ) 



4Q e 

3c 3 m £2 3c 2 m El 

a s = (cosd £2 )H (cosd E1 ) 

4C 



4C, 



dc 



dc 



3m E 2 

A'o = -(/i cos8 E 2 + f 2 sm8 E 2 ) 



4Q C 

3 m E,l 

4C. 



(f 3 cos8 El +f 4 sin8 E1 ) 



dc 



K w =-^(sm8 + co&8) , K n =-^- 
b l m E 2 d l m E j 

K \2 = ( C0S S E,2 + sin $E,2 ) + -7^- ( C0S S E,1 + Sin S E,\ ) 
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d,V, 



s ,e,\ = -T^icosS^ + sin<J £1 ), 

K s,E,2 =^ L ( cos ^,2+ sin ^,2) 



K 



d l m E ,V l 



dc 



s,SE,\ 



K 



s,SE,2 



K qd ~ 



2C se 
2C« 



(cos<5 £ j -sin<5 £ j) 
(cosc5 £ 2 -sin J £ 2 ) 



W £ 2 COS <J £ 2 



W £ j COS (5 £ [ 



f^y^- | m£ 2COS( J £ 2 + 



r d x +d 2 



2 

b ]+ b 



m E j cos 8 E j 



— |m £>2 sin£ £ 2 + 



m E,\ sin ^£,l 



V, 



b\+b 2 \ ( d 1 + d 2 . 

' 'm E2 smd E 2 +| — — — - \m F , sind F 



' l £,l sm ^£,1 



! ! w E 2 cos S E2 +\ 1+ ^ 2 |m £1 cos£ £] 



K„ 



2 



a: 



?,£,2 



2 

fc,+2>2 



a: 



a: 



q,SE,2 



2 

d l +d 2 



^ e COsJ £)2 



a: 



P,E,\ 



m E,\VdcS™8 El 



m E,2 V dc sin ^E,2 



''' ^ \KV,j-X' d I t Jco S S E , 



+ {V t , q +X q I t4 )^m5 EA W l 



dc 



K P , E ,2 = (V t , q + X q I t J^-)v dc sin«? 



£.2 



-(a 1+ a 2 )(V tJ -X d I tiq )V bS mS 



K 



p,SE,l ~ 



d t +d 2 



KV,, q + X q I t ,d)cosS EA 



(Vt,d-X' d I t JsinS Etl ]m Etl V dc 



K 



p,SE,2 



(v,, q + x q i tM ) 



r b x +b 2 



m E,2 V dc COS S E,2 



b l +b 2 



K, 



3d 3 m Ea V dc 



c,E,\ 



8C. 



(V,, d -X' d I t , q )m E 2 V dc smS E2 

(cos S E1 cos 8 E 2 + sin £ £ [ sin J £2 ) 



3^2™£,l V <fc 



8C 



(cos <5 £ j 2 + sin S E1 2 ) + w x 



dc 



3b 3 m E2 V d 2 . 2 , 

^,£,2= — ^ (cosc? £2 +sin<? £>2 ) 



8Q C 



3b 2 m E y dc 

H - (cos £ £>2 cos S E l + sin (5 £ 2 sin S El ) + w 2 



K 



c,SE,\ 



8Q C 
3</ 3 m £i2 m £il V dc 



8C. 



(-sin£ £J cos<J £ 2 + cos£ £1 sinc> £ 2 ) 



dc 



M 2 m EA V dc 



8C, 



(-sin£ £ j cosS E [ + cos <J £ j sin <J £ j ) + >v 3 



3*3'%,2 V<fc 



c,<5E,2 



8C. 



(- Sin <5 £ 2 cos <^£ 2 + cos $E 2 s i n <^£ 2 ) 



dc 



3b 2 m E1 m E2 V dc 



8C. 



(- sin £ £ 2 cos ^E,l + cos $E 2 s i n $E 1 ) + ^4 



a: 



v,£,l 



fifj +(i 2 



X q sin^ £>1 - 



d l +d 2 



X' d cosS EA ]V dc 



K, 



fV r,d)( b x +b 



v,E,2 



q" dc *> illu E,2 



K 



v,SE,\ 



d l +d 2 \ v > 



K 



v,SE,2 



2 



X' d m Etl V dc sinS EA + 



( 6> +b 2 




2 




%d 


+ 







X 'd V dc COsS E,2 



X q m 6 



\ +b 2 



q" l E,2 y dc ^ au E,2 



+ \ ^^\x' d m Ea V dc sinS Et2 



fl = | ~W \ m E,2 sin ^£,2 +1 \ m E,l Sin<? £,l 



h =^yJ' M £,2 COS<? £,2+^Y I"'/ 1 cos ^£.l 

f\ = I y j ™£,2 sin <?£,2 + 1 y | TO £,i sin <? £>1 



'_3 " 

^ 3 > 



7 £, W cos ^£,i + -7^ ^ sin <?£,i 



' 3m £ / 
^ 3m £t2 ^ 

v 4 Qc y 



7 £,2,d COsJ £i2 + 



7 £,1, 9 C0S ^£,1- 



7 £,2, 9 CO S ^£,2 



^ 3m £1 

v 4C dc J 



r ^E,2^ 



hxq Sill<? £,2 



hxd sin <?£,i 



v 4 Q C y 



7 £ ; 2,d sin<5 £>2 



APPENDIX C 



The investigated power system parameters used in 
simulation process are as follows: 

^£,1 = 30° ; <? £>2 = -25° ; X E l = X E 2 = 0.2 pu ; 
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X se = X te = O.lpu ; X bv =X d = 1.0 pu ; X q = 0.6pu ; 

X d ' = 03pu ; C dc = C se = 3.0/7M ; m E , = m £2 = 1 ; 
V Jc = = 2.0/w ; V b = \.0pu ; M = 2H = 8.0s ; 

/, = 50Hz ; T do ' = 5.044.5 ; T a = 0.05s ; K a = 50 ; D m = 
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